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Abstract 

Mid- Wave  Infrared  (MWIR)  semiconductors  are  of  use  to  the  Air  Force  for 
several  applications.  Ultrafast  spectroscopy  can  be  used  to  better  quantify  the  ef¬ 
fectiveness  of  semiconductor  opto-electronic  devices.  The  objective  of  this  thesis  was 
to  improve  the  procedure  for  taking  ultrafast,  time-resolved  measurements  of  photo¬ 
luminescence  from  MWIR  semiconductors.  Previous  work  has  used  a  mode-locked 
titanium  sapphire  (Ti:Saph)  laser  to  excite  the  semiconductor  sample  and  to  upcon- 
vert  the  photoluminescence  from  the  semiconductor.  Work  completed  in  this  thesis 
improved  on  the  techniques  developed  during  previous  work.  A  synchronous  Optical 
Parameter  Oscillator  (OPO)  will  be  used  to  convert  the  Ti:Saph  laser  (0.830  //m) 
into  1 . 3- /irn  signal  and  2. 3- /irn  idler  beams.  Whereas  the  Ti:Saph  reference  and  signal 
beams  were  too  energetic  to  excite  the  quantum  well  (QW)  layers  of  certain  MWIR 
opto-electronic  semiconductor  structures,  the  synchronous  OPO  allows  investigation 
of  the  active  layer  by  directly  exciting  the  QWs.  Photoluminescence  was  detected 
from  an  indium  arsenide  antimonide  (InAsSb)  sample  with  a  peak  wavelength  of  3.8 
/irn.  A  detailed  procedure  for  setup  of  the  TRPL  experiment  using  the  synchronous 
OPO  is  provided. 
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Ultrafast  spectroscopy  of  mid-infrared  semiconductors 


USING  THE  SIGNAL  AND  IDLER  BEAMS  OF  A  SYNCHRONOUS  OPTICAL 

PARAMETRIC  OSCILLATOR 

I.  Introduction 

With  the  introduction  of  the  FIM-92  Stinger  (Redeye  II)  by  General  Dynamics 
in  1967  [16],  there  has  been  a  persistant  threat  of  a  man-portable  air  defense  system 
(MANPADs)  attack  against  military  and  civilian  aircraft.  To  ensure  the  safety  of 
U.S.  civilians  and  military  personnel,  a  method  of  deterring,  detecting  and  defeating 
these  infrared  (IR)  based  systems  is  needed. 

Currently,  it  is  estimated  that  there  are  over  700,000  [11]  active  MANPADs  in 
the  world,  with  between  50,000  to  100,000  [4]  available  to  non-government  agencies. 
The  cost  of  a  single  unit  is  between  $35,000  and  $200,000  [9],  making  them  easily 
attainable  by  extremist  groups. 

As  a  result  of  the  threat  of  these  IR  systems,  work  has  been  done  to  reduce  the 
proliferation  of  these  threats  and  their  use  specifically  to  civilian  aircraft.  [18, 19,25] 
Two  methods  that  have  been  used  to  reduce  the  threat  have  been  purchasing  the 
outstanding  arms  to  deter  use  and  using  flares  to  defeat  them. 

Traditional  countermeasure  methods  against  IR  systems  have  been  employing 
Magnesium- Teflon- Viton  (MTV)  flares.  These  flares  work  by  deploying  a  pyrophoric 
material  that  ignites  in  contact  with  air.  Since  the  pyrophoric  material  is  hazardous, 
a  suitable  alternative  is  required.  Also,  MTV  flare  decoys  have  been  suitable  for 
first  generation  surface-to-air  missiles  (SAMs),  such  as  the  SA-7,  but  lack  the  ability 
to  counter  second  generation  SAMs.  [5]  First  generation  missiles  used  a  single  IR 
frequency  to  acquire  a  target.  This  acquisition  method  could  be  fooled  by  the  MTVs. 
Current  IR  threats  acquire  their  target  using  two  frequencies,  only  one  of  which  is 
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in  the  IR.  The  second  signal  is  detected  in  the  ultra-violet  (UV)  region.  MTV  flares 
have  a  low  success  rate  against  current  IR  missiles  because  of  the  UV  frequency  used 
to  supplement  target  acquisition. 

1.1  Motivation 

Two  current  technologies  are  being  investigated  to  defeat  IR  threats.  They  are 
the  U.S.  Army’s  Advanced  Threat  Infrared  Countermeasure  (ATIRCM)  and  the  U.S. 
Special  Operations  Command’s  Directional  Infrared  Countermeasure  (DIRCM).  [17] 
These  systems  employ  IR  lasers  broadcasting  several  jamming  algorithms  intended  to 
confuse  an  incoming  missile.  The  current  challenge  is  to  build  smaller  and  more  pow¬ 
erful  lasers.  One  such  program  underway  is  the  Defense  Advanced  Research  Projects 
Agency’s  (DARPA)  Steered  Agile  Beam  (STAB).  STAB’s  goal  is  to  develop  a  higher 
power  IR  laser  and  a  faster  moving  tracking  and  engagement  pod.  [13] 

Interband  IR  semiconductor  lasers  operate  using  the  decay  of  excited  electrons 
from  the  conduction  band  of  a  semiconductor  to  the  valence  band.  This  decay  results 
in  a  source  of  Mid- Wave  Infrared  (MWIR)  radiation.  Since  heat  signatures  from 
aircraft  naturally  broadcast  in  this  window,  MWIR  sources  are  ideal  for  defense  use. 
They  are  usable  in  an  atmospheric  electromagnetic  transmission  window  (2-5  /irn) . 
Figure  1.1  shows  the  transmittance  of  wavelengths  up  to  15  /mi  in  the  atmosphere. 
As  seen  in  Figure  1.1,  UV  lasers  are  impractical  as  a  countermeasure  due,  in  part,  to 
large  absorption  in  the  atmosphere  by  water  molecules. 

The  needs  for  our  aircraft  to  be  protected  from  hostile  munitions,  and  for  our 
munitions  to  be  assured  of  hitting  their  target,  means  the  development  of  IR  based 
counter  and  counter-counter  measures  are  interests  of  national  security.  It  is  these 
needs  that  are  driving  current  programs  investigating  and  developing  MWIR  opto- 
clectric  devices. 
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Figure  1.1.  Graph  of  atmospheric  transmittance  for  wavelengths  up  to  15  /irn  [1] 

1.2  Background 

A  better  understanding  of  MWIR  structures  will  aid  in  the  development  of 
higher  power  IR  lasers.  Three  areas  that  will  be  explored  a  bit  here  and  in  greater 
depth  in  Chapter  2  are  the  rate  equation,  photoluminescence  (PL)  and  time-resolved 
photoluminescence  (TRPL).  The  rate  equation  describes  the  decay  of  electrons  from 
the  conduction  band  to  the  valence  band  within  a  MWIR  semiconductor  structure. 
PL  is  the  term  for  the  light  emitted  in  this  decay  when  the  electrons  are  excited  by 
optical  radiation.  TRPL  is  the  process  that  will  be  used  to  study  PL. 


1.2.1  Rate  Equation.  For  bulk  semiconductors,  carrier  dynamics  can  be 
described  by  the  rate  equation 

1  AR 

(1.1) 


1 

r 


5R 

5n 


where  r  is  the  differential  carrier  lifetime,  R  is  the  recombination  rate  and  n  is  the 
excited  carrier  density.  A  power  series  solution  can  be  assumed  for  R  of  the  form 


R{n )  =  Ain1. 


(1.2) 


i= 0 
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Substituting  Equation  (1.2)  into  (1.1)  yields 


1  OO 

-  =  J2Anl~1-  (1-3) 

2=1 

An  expansion  of  Equation  (1.3)  would  involve  an  infinite  number  of  terms  and 
would  be  impossible  to  fit  to  for  any  computer.  As  a  result,  the  generally  accepted 
method  of  dealing  with  this  problem  is  to  truncate  all  but  the  first  three  terms  and 
to  assume  all  other  terms  are  negligible  in  comparison.  This  assumption  leads  to 

—  ^  =  A  ■  n  +  B  ■  n2  +  C  ■  n3.  (1.4) 

This  is  the  rate  equation  for  recombination  within  a  semiconductor.  The  coefficients 
A,  B  and  C  are  related  to  the  probability  of  relaxation. 

1.2.2  Photoluminescence.  PL  is  light  emitted  by  a  sample  following  the 
absorption  of  optical  photons.  It  occurs  when  an  electron  is  excited  to  a  higher 
energy  level  by  absorbing  a  photon,  and  then  spontaneously  decays  to  a  lower  energy 
level,  emitting  a  second  photon  in  the  process.  To  conserve  energy,  the  emitted  photon 
cannot  have  more  energy  than  the  exciting  photon.  [21] 

1.2.3  Time  Resolved  Photoluminescence.  TRPL  is  a  procedure  that  will  be 
used  to  study  the  PL  produced  by  a  semiconductor.  It  involves  a  second  non-collinear 
laser  beam  to  sample  the  PL.  Using  a  nonlinear  crystal,  the  sample  beam  can  combine 
with  the  PL  through  nonlinear  effects;  this  combination  is  called  upconversion.  Figure 
1.2  shows  how  a  sample  beam  combines  with  the  PL  to  yield  a  third  beam. 

When  the  sampling  beam  is  pulsed,  time  resolution  is  possible.  Using  a  pulsed 
laser,  the  PL  can  be  sampled  at  varying  times.  Figure  1.3  shows  how  this  is  accom¬ 
plished.  A  pulse  at  t=0  excites  the  semiconductor  which  then  begins  to  luminesce. 
The  intensity  can  then  be  resolved  temporally  by  delaying  the  sampling  pulse.  To 
determine  the  time  at  which  t=0,  it  is  necessary  to  find  zero  path  length  (ZPL).  ZPL 
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is  defined  as  the  path  length  difference  between  the  sampling  pulse  and  the  excitation 
pulse/PL  such  that  the  start  of  the  excitation  pulse/PL  overlaps  with  the  sampling 
pulse  in  the  nonlinear  crystal. 


Figure  1.2.  A  diagram  of  upconversion.  Two  input  beams  enter  the  nonlinear  crys¬ 
tal  from  the  right  and  result  in  a  third  upconverted  beam  on  the  left.  The  beams 
above  and  below  the  upconverted  beam  are  the  input  beams  if  there  is  no  upcon¬ 
version.  Due  to  conservation  of  momentum  the  resultant  beam  will  be  able  to  be 
separated  from  the  PL  and  sample  beam  components. 


Excitation  Pulse 

j  Luminescence 

*  t  Sampling  Pulse 


Variable  Delay 


Figure  1.3.  A  diagram  of  TRPL.  A  pulse  at  t=0  excites  the  semicondutor,  which 
then  begins  to  luminesce;  this  is  shown  as  the  luminescence.  The  sampling  pulse  can 
then  be  delayed  and  upconverted  with  the  PL  to  obtain  time  resolution  of  the  PL 
intensity. 
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1.3  Previous  Work 


There  have  been  several  attempts  to  achieve  TRPL  at  the  Air  Force  Institute 
of  Technology  (AFIT)  beginning  in  1997  with  Cooley.  [7]  Each  experiment  has  added 
to  the  understanding  of  TRPL.  Cooley  began  by  investigating  the  Shockley-Read- 
Hall  (SRH),  radiative  and  Auger  recombination  rates  of  indium  aluminium  arsenide 
antimonide  (InAlAsSb)  multiple  quantum  well  (MQW)  diode  laser  structures  emitting 
near  3.3  /jrn,  and  showed  that  the  SRH  rate  dominated  recombination  at  temperatures 
between  77K  and  150K.  Gorski  [10]  demonstrated  the  Auger  rates  were  as  low  as 
reported  by  Cooley,  and  showed  that  increases  in  the  InAs/GalnSb  well  width  of 
a  type-II  MQW  decreased  the  recombination  time.  Cumblidge  and  Johnson  [8, 12] 
showed  a  curve  fit  of  the  PL  data  that  forced  one  of  the  three  terms  in  the  rate 
equation  to  be  non-physical. 

This  dissertation  and  these  theses  used  a  mode- locked  titanium  sapphire  (Ti:Saph) 
laser  operating  at  ~  800  nm.  The  linearly  polarized  light  was  rotated  to  45°  and  split 
by  a  polarizing  beam  splitter  into  horizontally  and  vertically  polarized  light.  The 
excitation  of  the  semiconductor  and  the  upconversion  was  provided  by  this  800-nm 
light.  In  the  case  when  the  cladding  would  have  absorbed  the  Ti:Saph  laser,  the  clad 
was  removed  to  observe  the  carrier  dynamics  of  the  MQWs.  The  drawback  of  these 
studies  was  that  the  full  device  could  not  be  studied.  Roberts  [20]  introduced  an 
optical  parametric  oscillator  (OPO)  to  attempt  to  excite  the  core  of  type-II  MQWs 
without  removing  the  clad;  however,  he  was  only  able  to  achieve  time-averaged  PL 
with  a  Spectra-Physics  Opal  OPO. 

1.4  Approach 

In  an  interband  IR  semiconductor,  recombination  after  an  impulse  excitation 
does  not  occur  instantaneously,  rather  it  occurs  over  a  period  of  nanoseconds  [10] .  To 
develop  lasers  capable  of  delivering  additional  power  on  target,  the  electrical  proper¬ 
ties  of  semiconductors  need  to  be  investigated.  Studied  semiconductors  include  intrin- 
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sic  InAs,  AlAso.085Sbo.915/I1io.88Alo.12Aso.878Sbo.122  tvpe-I  MQW  and  InAs/Gao.ssbio.isSb 
type-II  MQW.  [8] 

As  stated  in  Section  1.3,  previous  work  used  a  mode-locked  Ti:Saph  laser  operat¬ 
ing  at  0.800  /iin  and  removed  the  clad  in  the  semiconductor  opto-electronic  structure. 
The  procedure  used  in  this  thesis  retains  the  clad  and  adds  a  synchronous  OPO  which 
allows  0.832-/mi  light  from  the  Ti:Saph  to  be  converted  to  1.300-/mi  and  2.311-yum 
beams.  The  2.311-/mi  beam  is  used  to  excite  the  core  of  the  semiconductor  structure 
and  the  1.3()0-/jrn  is  used  as  the  sampling  pulse. 

For  this  procedure  to  work,  a  technique  to  set  up  and  align  the  equipment 
needs  to  be  established.  A  thermal  camera  and  an  IR  sensitive  card  will  aid  in  the 
establishment  of  ZPL.  Following  an  attempt  at  ZPL,  PL  was  detected  from  a  MWIR 
semiconductor. 

1.5  Problem  Statement 

Since  laser  development  is  an  exercise  in  efficiency,  the  study  of  carrier  dynamics 
of  potential  semiconductors  is  an  important  step  in  determining  effectiveness  of  a 
semiconductor  laser.  Observing  the  temporal  decay  of  the  PL  of  an  excited  sample 
will  allow  the  temporal  decay  of  the  carrier  densities  to  be  determined.  From  these 
densities,  the  coefficients  of  the  rate  equation  (A,  B  and  C)  can  be  fit. 

Previous  research  investigated  carrier  dynamics  with  a  higher  photon  energy 
Ti:Saph  laser.  This  study  will  introduce  an  OPO  to  reduce  the  photon  energy  of  the 
excitation  pulse  to  better  examine  the  QW  carrier  dynamics.  R  will  also  produce 
documentation  and  a  procedure  for  the  setup  of  upconversion  to  be  used  for  future 
study  of  TRPL. 

1.6  Results 

This  thesis  provides  a  step-by-step  process  for  aligning,  using  and  calibrating 
a  Ti:Saph  laser  and  a  synchronous  OPO;  establishing  sampling,  excitation  and  PL 
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pathlengths;  collection  of  PL;  selection  of  a  nonlinear  crystal  for  upconversion;  and 
attempting  to  find  ZPL  for  the  collection  of  the  npconverted  signal  for  analysis. 

1.7  Overview 

This  thesis  is  broken  up  into  five  chapters  and  one  appendix.  Chapter  Two 
discusses  the  theory  of  carrier  dynamics.  Chapter  Three  describes  the  experiment. 
Chapter  Four  details  the  results  and  Chapter  Five  is  the  conclusion.  The  appendix 
includes  a  discussion  of  software  and  the  Mathcmatica®  code  used  in  this  thesis. 


II.  Theoretical  Formulation 


The  lasing  efficiency  of  a  MWIR  semiconductor  structure  is  highly  dependent  on 
carrier  dynamics.  Besides  radiative  decay,  there  are  also  non-radiative  effects  that 
heat  the  material  and  reduce  efficiency.  A  better  IR  laser,  using  a  MWIR  material, 
requires  better  characterization  of  these  non-radiative  effects.  Understanding  these 
effects  allows  semiconductor  crystal  growers  to  design  better  structures  to  mitigate 
the  non-radiative  transitions.  This  would  lead  to  higher  power  IR  lasers  suitable  for 
defense  applications. 

2.1  Background 

2.1.1  Quantum  Wells.  MWIR  structures  are  manufactured  by  epitaxial 
growth  on  a  crystalline  structure.  There  are  three  kinds  of  epitaxial  structures  for 
semiconductors:  types-I,  -II  and  -III  MQWs.  Type-I  MQWs  are  characterized  by 
having  the  QW  in  the  conduction  and  valence  bands  completely  contained  in  a  single 
layer.  Type-II  MQWs  are  created  when  the  QW  in  the  conduction  band  is  in  one 
layer  and  the  QW  in  the  valence  band  is  in  the  adjacent  layer.  Finally,  a  type-III 
MQW  is  created  when  the  conduction  band  edge  of  one  layer  has  a  lower  energy 
than  the  valence  band  edge  of  the  adjacent  layer.  Carriers  can  still  transition  from 
the  conduction  band  to  valence  band  as  the  allowable  conduction  band  energies  are 
higher  than  the  allowable  valence  band  energies.  Figures  2.1  (a)-(c)  illustrate  the 
energy  gaps  as  a  function  of  location  in  the  structure. 

2.1.2  Electron- Hole  Pair  Formation  and  Recombination.  Photon  absorption 
works  by  taking  the  energy  from  a  photon  to  promote  an  electron  over  the  band  gap. 
This  process  is  instantaneous.  This  absorption  creates  an  electron-hole  pair,  which 
is  generally  excited  deep  into  the  bands.  The  electron  and  hole  then  relax  through 
carrier  cooling  to  the  band  edge.  This  process  is  on  the  order  of  picoseconds.  Finally, 
the  electron-hole  pair  relax  in  the  form  of  recombination.  This  last  process  occurs 
over  a  period  of  nanoseconds. 
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Figure  2.1.  Energy  diagram  for  types-I  (a),  -II  (b)  and  -III  (c)  MQW.  The  outer, 
protective  layer  is  known  as  the  cladding. 

Carrier  recombination  is  characterized  by  three  types  of  transitions:  radiative, 
trap  or  defect,  and  Auger.  Defect  is  one  of  the  undesirable  transitions.  Figure  2.2(a) 
illustrates  how  an  excited  electron  can  decay  from  the  conduction  band  to  the  valence 
band  in  several  steps.  Each  step  is  caused  by  impurities  in  the  semiconductor  resulting 
in  energy  levels  within  the  designed  band  gap.  This  releases  several  phonons  instead 
of  the  desired  single  photon.  Radiative  recombination  is  the  preferred  transition  as 
it  releases  a  photon  of  approximately  the  band  gap.  When  the  sample  is  excited  by 
optical  radiation,  this  is  PL.  Figure  2.2(b)  shows  this  recombination.  Finally,  Auger 
recombination  is  defined  as  a  three-body  transition.  It  involves  an  electron  in  the 
conduction  band  transitioning  to  the  valence  band.  The  energy  from  this  transition 
is  not  released  as  a  photon,  but  rather  is  used  to  excite  a  second  electron  from  deep  in 
the  valence  band  to  another  valence  band  state.  This  second,  now  excited,  electron  is 
then  free  to  return  to  thermal  equilibrium  without  releasing  a  photon.  This  process 
is  shown  in  Figure  2.2(c).  This  is  the  second  undesirable  form  of  recombination. 

2.2  Theory  of  carrier  dynamics 

Once  an  electron-hole  pair  is  generated,  the  system  is  no  longer  in  thermal 
equilibrium.  However,  as  in  section  2.1.2,  the  intraband  relaxation  of  the  electrons 
and  holes  occurs  over  picoseconds;  the  conduction  band  electrons  can  be  considered 
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Figure  2.2.  Illustration  of  the  three  forms  of  recombination:  defect  (a),  radiative 
(b)  and  Auger  (c).  Radiative  recombination  is  the  preferred  recombination. 


in  equilibrium  among  themselves,  as  are  the  valence  band  holes.  This  situation  is 
known  as  quasi-equilibrium. 

The  development  for  calculating  the  carrier  density  begins  with  the  density  of 
states  for  the  electrons  and  holes,  which  are  defined  as 


QC(E )  = 

(2m*)3/2  , _ 

Qv(e)  =  VE^EdE. 


(2.1) 


where  h  is  the  Plank  constant  divided  by  27r,  m^2  and  are  the  effective  masses 

of  the  electrons  and  holes,  respectively,  E  is  the  energy,  and  Ec  and  Ev  are  the  band 
edge  energies  for  the  conduction  and  valence  bands  respectively.  The  first  equation  is 
the  density  of  states  for  the  electrons  in  the  conduction  band  and  the  second  is  the 
density  of  states  for  the  holes  in  the  valence  band.  [15]  Fermi-Dirac  statistics  provide 
the  next  step  as  the  distribution,  fn(E),  represents  the  probability  a  conduction-band 
energy  level,  E ,  is  occupied  by  an  electron.  fp(E)  then  represents  the  probability 
a  valence-band  energy  level,  A,  is  occupied  by  a  hole.  Equation  (2.2)  shows  the 
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respective  distribution  functions. 


i  +  e(Ec-Efc)/kT 
fP(E)  =  1  -  1  +  e(E}v-Ev)/kT 


(2.2) 


Efc  and  Efv  are  the  quasi-Fermi  levels  for  the  conduction  and  valence  bands,  respec¬ 
tively. 

The  carrier  densities  are  then  n(E)  =  gc(E)  ■  fn(E )  for  electrons  and  p(E)  = 
Qp(E)-fp(E)  for  holes.  [15,24]  Figure  2.3  illustrates  a  semiconductor  in  quasi-equilibrium. 
It  shows  the  two  distribution  functions  and  the  resulting  carrier  densities  in  each  band. 
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Figure  2.3.  A  semiconductor  in  quasi-equilibrium.  fn{E)  is  the  probability  that  a 
conduction-band  energy  level  E  is  occupied  by  an  electron.  The  probability  that  a 
valence-band  energy  is  occupied  by  a  hole  is  1  -  fp(E).  n(E)  and  p(E)  are  obtained 
by  multiplying  the  Fermi  function  by  the  density  of  states.  [21] 


2.2.1  Carrier  densities  from  intensity.  The  rate  equation  now  needs  to  be 


related  to  a  measurable  value.  Since  the  PL  signal  can  be  measured  and  not  the 
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carrier  density  directly,  the  rate  equation  needs  to  be  related  to  the  PL.  From 


(2.3) 


the  relation  between  PL  intensity  and  carrier  density  can  be  established.  [7]  n(io)  is 
the  initial,  and  maximum,  carrier  density  in  the  QW  at  excitation,  e  is  the  energy 
per  pulse  and  is  divided  by  the  energy  per  photon,  h  ■  v.  Nqw  is  the  number  of  QWs. 
wp  and  wpi  are  the  waist  of  the  pump  and  PL,  respectively.  R  is  the  Fresnel  reflection 
at  the  semiconductor  surface,  a  is  the  absorption  coefficient.  Finally,  ts  and  tsAR 
are  the  thicknesses  of  the  QW  and  barrier,  respectively,  within  the  semiconductor 
structure.  Equation  (2.3)  is  used  to  calculate  the  initial  carrier  densities  at  various 
excitation  powers  and  is  plotted  against  the  experimentally  measured  PL  intensities 
at  the  same  excitation  powers.  This  curve  is  shown  in  Figure  2.4.  [10] 
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Figure  2.4.  A  plot  of  carrier  density  as  a  function  of  luminescence.  [10] 


The  equation  used  to  fit  the  data  in  Figure  2.4  is  given  by 


n(L)  =  u\  ■  L  +  U2  ■  LU3 


(2.4) 
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where  U\,  w2  and  w3  are  coefficients  that  provide  the  best  fit.  Using  TRPL,  the  carrier 
densities  can  then  be  plotted  as  a  function  of  time. 

2.3  Theory  of  Upconversion 

TRPL  is  contingent  on  obtaining  upconversion  of  the  PL.  Upconversion  is  the 
process  of  sum  frequency  generation  (SFG).  SFG  occurs  when  two  optical  waves  of 
different  frequencies  are  mixed  to  generate  a  third  wave  of  a  higher  frequency.  This 
mixing  occurs  when  the  energy  conservation  and  phase  matching  conditions  are  met. 
The  energy  conservation  condition  is 


Ui+U2  =  u3.  (2.5) 

where  uqj2  are  the  frequencies  of  the  input  beams  and  cu3  is  the  frequency  of  the  output 
beam.  The  phase  matching  condition  is 

ki  +  k2  =  k3  (2-6) 

where  k12,3  are  the  propagation  vectors  of  the  input  and  output  beams,  respectively. 
The  difficulty  obtaining  SFG  comes  from  simultaneously  solving  Equation  (2.5)  and 
(2.6). 

SFG  will  not  occur  in  an  isotropic  material;  it  will  only  occur  in  anisotropic  me¬ 
dia.  There  are  two  kinds  of  anisotropic  crystals:  uniaxial  and  biaxial.  Biaxial  crystals 
have  three  different  principle  indices  of  refraction.  Due  to  symmetry  of  the  crystalline 
structure,  uniaxial  crystals  only  have  two  indices  of  refraction  known  as  an  ordinary 
index,  n0,  and  an  extraordinary  index,  ne.  A  positive  uniaxial  crystal  has  ne  >  n0 ; 
a  negative  uniaxial  crystal  has  ne  <  na.  The  different  indices  of  refraction  raise  an 
additional  challenge  to  phase  matching  as  polarization  also  needs  to  be  considered. 
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Using  a  negative  uniaxial  crystal  as  an  example,  the  phase  matching  condition 


for  the  collinear  case  becomes 


ne{uj i)  •  u>i  +  na( u2)  ■  u2  =  ne(u3 )  •  uj3  (2.7) 

An  example  of  a  negative  uniaxial  crystal  is  lithium  iodate  (LilCU).  This  crystal  is 
important  because  it  allows  phase  matching  for  ZPL.  As  Figure  1.2  showed,  ki  and 
k2  are  not  collinear.  This  results  in  ne  being  a  function  of  incident  angle,  as  well 
as  uj.  However,  the  complications  of  SFG  are  simplified  by  the  use  of  two  software 
programs:  SNLO  nonlinear  optics  code  and  a  Mathematica  notebook.  Using  these 
two  programs,  the  numerical  challenges  are  alleviated.  Additional  discussion  of  these 
programs  is  included  in  Appendix  A. 

2.4  Summary 

Carrier  dynamics  and  SFG  were  discussed  in  this  chapter.  TRPL  uses  SFG  to 
provide  time  resolution  of  PL  decay.  Relating  the  carrier  densities  to  the  intensity  of 
the  PL  provides  time  resolved  carrier  dynamics.  The  resultant  data  can  then  be  fit 
to  the  rate  equation  to  obtain  the  recombination  coefficients. 
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III.  Experiment 

This  chapter  discusses  the  experimental  setup,  equipment  used,  and  the  documented 
procedure  for  TRPL. 

3.1  Experimental  Procedure 

TRPL  spectroscopy  was  first  described  by  Mahr  [14]  and  detailed  by  Shah.  [22] 
It  allows  for  the  characterization  of  PL  on  a  time  scale  of  the  same  order  as  the 
excitation  pulse.  The  intensity  of  the  upconverted  PL  pulse  is  related  to  the  carrier 
densities  via  the  relationships  as  discussed  in  Chapter  2.  These  carrier  densities  are 
then  used  to  fit  theoretical  models  describing  radiative  and  non-radiative  relaxation 
processes. 

3.2  Experimental  Setup 

This  thesis  uses  a  mode-locked  Ti:Saph  laser  and  a  synchronous  OPO  to  generate 
the  photons  of  proper  wavelength  necessary  to  excite  the  core  of  various  MWIR  semi¬ 
conductor  structures.  There  are  four  distinct  optical  paths:  sampling  (signal  beam), 
excitation  (idler  beam),  luminescence  (PL),  and  upconverted.  Figure  3.1  shows  the 
setup  of  this  experiment.  The  sampling  path  has  a  fixed  optical  path  length  (OPL) 
and  has  a  focus  at  the  non-linear  crystal.  The  second  path,  the  excitation  path,  trav¬ 
els  a  variable  OPL  and  is  incident  on  a  liquid  nitrogen  cooled  semiconductor.  The 
variation  of  this  path  is  provided  by  a  translatable  mirror  which  is  able  to  temporally 
shift  the  luminescence  relative  to  the  arrival  of  the  sampling  pulse.  By  shifting  the 
start  of  PL,  the  sampling  path  is  able  to  experience  upconversion  with  the  lumines¬ 
cence  path  -  the  path  from  the  semiconductor  -  at  the  non-linear  crystal  at  various 
times  along  the  PL  lifetime.  The  fourth  path  is  the  upconverted  signal  as  it  exits  the 
non-linear  crystal  and  is  directed  to  a  photomultiplier  tube  (PMT). 

3.2.1  Pump  laser.  The  first  component  of  this  setup  is  a  Coherent  Verdi 
Nd:YVC>4  laser  operating  at  a  wavelength  of  532  nm.  The  operating  power  is  12.5 
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Figure  3.1.  Setup  of  the  TRPL  equipment 
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W.  It  uses  a  water  cooler  as  a  heat  sink.  To  prevent  damage  to  the  laser,  the  cooler 
needs  to  be  at  18°C;  cooling  can  take  up  to  15  minutes  to  reach  steady  state. 

3.2.2  Ti:Saph  laser.  [6]  The  Verdi  enters  a  Coherent  Mira  Ti:Saph  femtosec¬ 
ond  laser.  This  laser  amplifies  wavelengths  between  680  and  1100  nm  and  further 
restricts  the  wavelengths  to  150-nm  ranges  centered  on  three  wavelengths.  Each  las¬ 
ing  wavelength  requires  that  the  cavity  be  an  integer  number  of  wavelengths;  this 
is  known  as  a  longitudinal  mode.  It  utilizes  mode-locking  to  create  a  pulsed  laser 
operating  at  76  MHz  at  1.4  W.  Mode  locking  occurs  when  the  frequencies  generated 
undergo  destructive  interference  from  all  available  modes  excepting  at  the  pulse  where 
the  modes  constructively  interfere.  More  modes  available  yield  a  smaller  pulse  width 
so  a  high  modal  volume  is  encouraged.  The  Ti:Saph  crystal  absorbs  the  500-nm  light 
from  the  pump  and  reemits  is  at  the  available  frequencies. 

This  study  requires  the  OPO  to  receive  light  of  832  nm.  Since  there  is  a  water 
absorption  line  at  820  nm,  the  humidity  is  monitored  and  controlled.  There  is  a 
relative  humidity  monitor  included  on  the  controller  for  the  Ti:Saph  which  should 
be  below  25%  for  mode  locking  to  occur.  If  the  humidity  is  above  25%,  then  the 
humidity  can  be  controlled  via  a  nitrogen  tank  connected  to  the  Ti:Saph. 

For  ease  of  use  of  the  OPO,  maximum  power  is  required.  After  turning  on  the 
Verdi,  there  is  a  start-up  period  of  about  2  minutes  as  the  power  slowly  ramps  up  to 
the  operating  condition.  After  this  start  up  period,  the  power  for  a  given  wavelength 
should  be  maximized.  This  is  accomplished  by  tuning  the  pump  optic  mirrors,  the 
flat  end  mirror  and  the  slits,  all  as  described  in  the  manual. 

By  adjusting  the  bircfringent  filter,  the  wavelength  can  be  adjusted  between  680 
and  1100  nm.  Rotation  of  the  control  clockwise  decreases  the  wavelength  while  a  coun¬ 
terclockwise  rotation  increases  it.  As  the  wavelength  is  adjusted  in  small  increments, 
the  laser  should  be  kept  mode  locked  and  at  maximum  power.  The  wavelengths  that 
are  only  achievable  as  a  continuous  wave  should  be  avoided  for  use  in  the  OPO. 
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3.2.3  OPO.  [3]  The  830-nm  mode-locked  Ti:Saph  laser  enters  the  Coherent 
OPO  from  a  pair  of  mirrors  as  depicted  in  Figure  3.1.  Upon  entering  the  device,  the 
light  enters  a  ring  mirror  configuration.  Within  the  OPL  of  the  PP830-Ring  cavity, 
there  is  a  nonlinear  crystal.  When  under  a  stress,  it  downconverts  the  incoming 
Ti:Saph  laser  into  two  parts.  The  idler  and  signal  do  not  exit  the  OPO  at  the  same 
time;  the  signal  leads  the  idler  by  3  ns.  This  will  be  discussed  further  in  the  alignment 
procedure. 


3.2.3. 1  Discussion.  The  downconversion  of  the  Ti:Saph  is  an  impor¬ 
tant  process  to  understand,  as  the  idler  is  key  to  PL.  Upconversion  and  downconver¬ 
sion  require  the  conservations  of  energy  and  momentum.  Since  the  pump  Ti:Saph  is 
lasing  at  832  nm  and  the  OPO  signal  will  nominally  be  set  to  1300  nm,  from  Equation 
(3.1),  the  idler  will  be  at  2311  nm.  Since  the  idler  will  be  incident  on  the  semiconduc¬ 
tor,  the  photon  energy  needs  to  be  such  that  the  core  of  the  MWIR  device  structure 
will  be  excited  while  the  clad  remains  transparent  to  this  energy. 


3. 2. 3. 2  Operation.  The  OPO  outputs  three  lasers:  a  tunable  signal, 
a  feedback  signal,  and  an  idler.  To  begin,  the  feedback  signal  must  be  sent  into  a 
spectrometer;  however,  directing  the  full  intensity  of  the  laser  into  the  spectrometer 
will  irreparably  damage  it.  Instead,  a  beam  splitter  must  be  used  to  reduce  the 
irradiance.  Also,  due  to  the  path  in  the  OPO,  there  are  three  beams  exiting;  the 
strongest  of  the  three  is  to  be  directed  into  the  spectrometer.  This  feedback  provides 
the  OPO  software  a  response  to  the  adjustments  to  the  OPO  controller.  The  pump 
and  signal  will  be  detected  and  displayed  on  the  control  screen  when  the  path  is 
aligned.  Alignment  is  done  by  directing  the  strongest  beam  into  the  center  of  the 
focusing  lens  and  then  into  the  entrance  slit.  If  the  beam  does  not  enter  the  center 
of  the  focusing  lens  and  entrance  slit,  there  will  be  no  signal  displayed. 
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Since  tuning  for  the  bulk  crystal  is  conducted  via  pump  wavelength  tuning,  prior 
knowledge  of  the  Ti:Saph  laser  wavelength  is  needed  or  guess  work  will  be  necessary. 
The  piezo  voltage  applied  to  the  non-linear  crystal  can  be  modulated  between  0  and 
4095  in  arbitrary  units.  Based  on  the  incoming  wavelength,  the  voltage  needs  to  be 
approximately  tuned  based  on  Table  3.1.  In  addition  to  changing  the  piezo  voltage, 
for  A  signal  between  1100  and  1200  nm,  the  6- mm  dispersion  block  should  be  used. 
For  all  X signal  greater  than  1200  nm,  the  12-mm  dispersion  block  is  required.  Below 
A  signal  of  1100  nm,  no  dispersion  block  should  be  used. 


Table  3.1.  Variation  of  pump  wavelength  to  OPO  output  based  on  Mira  documen¬ 
tation  _ 


A  Signal  (nm) 

Power  (mW) 

A Ti-.SAFF  (nm) 

Piezo  Voltage  (A.U.) 

1100 

500 

840 

148 

1150 

560 

829 

500 

1200 

500 

829 

1100 

1250 

450 

829 

1300 

400 

830 

1700 

1350 

380 

832 

1748 

1400 

250 

831.5 

1450 

240 

831.5 

2148 

1500 

320 

831.5 

1550 

280 

831.5 

1600 

150 

836 

2678 

Variation  in  day-to-day  operation  of  the  equipment  requires  a  daily  tuning  of 
the  piezo  voltage  to  yield  the  desired  output  signal  wavelength.  The  controller  can 
automatically  control  the  signal  within  a  user-specified  tolerance  (up  to  AA  =  1  nm) 
provided  the  signal  from  the  feedback  is  strong  enough  to  be  detected.  If  the  feedback 
wavelength  is  weak,  the  input  from  the  Ti:Saph  needs  to  be  adjusted  as  detailed  above. 

The  controller  is  an  APE  OPO  PP  Automatic  supplied  from  Coherent  with  the 
OPO.  The  different  display  options  are  shown  under  “Spectrum  Menu”.  From  here, 
the  pump,  signal,  and  a  “vis”  wavelength  can  be  displayed.  Specific  to  this  thesis, 
the  “vis”  was  not  seen  and  was  not  needed.  [3] 
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The  signal  can  be  tuned  from  1100  -  1600  nm  based  on  the  input  pump  wave¬ 
length  and  piezo  voltage.  Setting  the  signal  is  done  by  slowly  adjusting  the  pump 
wavelength  while  keeping  a  mode-locked  laser  until  the  desired  pump  wavelength  is 
reached.  Once  the  Ti:Saph  is  set  to  the  proper  wavelength,  Table  3.1  gives  a  range  for 
the  piezo  voltage.  Caution  is  advised  when  operating  on  a  pump  wavelength  that  is 
barely  detectable,  as  fluctuations  in  operating  conditions  can  reduce  the  signal  below 
the  detection  threshold.  This  causes  the  automatic  piezo  control  to  attempt  to  find 
the  signal  again  by  increasing  the  piezo  voltage.  This  can  ruin  a  data  run  if  the  signal 
is  not  optimized.  The  signal  is  used  as  the  sampling  light. 

The  OPO  that  is  used  has  had  the  idler  option  installed.  The  idler  is  the  excess 
from  the  downconversion  process,  and  while  not  directly  monitored,  is  related  to  the 
input  and  signal  wavelengths  via  Equation  (3.1).  The  idler  is  used  to  excite  the 
semiconductor  structure’s  inner  active  region.  The  idler  beam  also  has  components 
from  the  pump  at  832  nm  in  addition  to  the  2311-nm  light.  The  832-nrn  component 
was  reduced  by  placing  a  long  pass  filter  in  the  path  of  the  idler. 

3.2.4  Alignment.  Alignment  for  this  procedure  is  crucial  to  success.  The 
OPLs  for  the  three  paths  are  on  the  order  of  meters.  As  a  result,  if  the  laser  is  not 
centered  on  the  center  of  the  mirrors,  alignment  becomes  challenging. 

3.2.4- 1  Zero  Path  length.  Since  the  signal  and  the  idler  have  different 
OPLs  inside  the  OPO,  they  exit  at  different  times.  Since  the  upconversion  process 
will  require  phase  matching,  the  OPLs  on  the  optical  table  must  be  established  such 
that  the  overall  OPL  of  each  path  is  the  same.  The  idler  leads  the  signal  by  10  ns  or 
the  signal  leads  the  next  idler  pulse  by  about  3  ns.  Based  on  the  configuration  of  the 
table,  the  idler  path  is  5.35  nr  at  a  minimum;  the  signal  path  is  2.53  m  at  a  minimum. 
The  signal  path  needs  to  be  either  0.9  m  longer  than  the  idler  path  or  2.998  m  shorter. 
Since  the  signal  path  is  2.82  m  shorter,  the  idler  path  can  be  increased  by  0.178  m  to 
get  ZPL.  Fine  adjustments  in  the  OPL  is  accomplished  with  the  translatable  stage. 
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Previous  work  used  an  IR  sensitive  card  from  Spectra-Physics  or  night  vision 
goggles  to  directly  observe  the  Ti:Saph  laser.  Since  both  the  signal  and  idler  have 
a  longer  wavelength,  neither  technique  worked  well  for  this  research.  The  solution 
to  this  problem  was  to  use  an  Inframetrics  760  thermal  viewer  and  a  Thor  labs  VC- 
VIS/IR  card.  The  IR  viewer  requires  approximately  a  10  minute  startup  time  to  allow 
the  camera  to  cool  down.  To  use  the  thermal  camera,  the  idler  was  allowed  to  heat  a 
piece  of  cardstock  to  a  temperature  detectable  by  the  IR  viewer.  The  challenge  of  this 
technique  is  that  the  response  time  is  approximately  20-30  seconds  of  heating/ cooling 
to  detect  variations  in  the  alignment.  This  technique  worked  well  when  the  mirrors 
were  across  the  table  from  each  other.  When  the  mirrors  were  on  the  same  side  of  the 
table,  the  VC-VIS/IR  card  with  the  night  vision  goggles  was  able  to  detect  the  signal 
only  in  a  darkened  room.  The  VC-VIS/IR  card  can  detect  the  signal  component 
directly.  The  difficulty  with  this  method  of  detection  is  that  if  the  power  from  the 
OPO  is  not  maximized,  by  having  the  power  from  the  Ti:Saph  maximized,  then  the 
intensity  reradiating  from  the  IR  card  will  be  weak. 

For  ZPL,  the  laser  pulses  must  be  phase  matched  at  the  LHO3  crystal.  In 
addition,  they  need  to  be  overlapping.  By  replacing  the  LiI03  crystal  with  a  pinhole, 
it  can  be  assured  the  beams  occupy  the  same  space  by  passing  the  idler  and  the  signal 
beams  through  the  pinhole.  A  photodetector  is  employed  for  alignment  on  the  other 
side  of  the  pinhole.  By  adjusting  the  focusing  mirror  directly  preceding  the  pinhole, 
and  passing  the  lasers  through,  the  intensity  in  the  photodetector  can  be  maximized. 
This  corresponds  to  a  focus  at  the  pinhole.  At  this  point,  the  paths  are  aligned  up  to 
the  non-linear  crystal. 

To  align  the  PL  path,  a  Helium  Neon  (HeNe)  laser  was  placed  in  line  with  the 
Triax  320  monochromator  as  in  Figure  3.1.  It  was  projected  between  the  off-axis 
parabolic  mirrors  towards  the  pinhole.  Next,  it  was  focused  through  a  63-mm  lens, 
through  the  pinhole,  where  the  signal  was  maximized  with  the  photodiode.  It  was 
then  projected  between  the  other  pair  of  off-axis  parabolic  mirrors  onto  the  cold  finger 
in  the  vacuum  chamber. 
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3.2.5  Procedures.  The  following  section  describes  the  methodology  to  turn 
on  the  equipment  used,  adjust  it  and  shut  it  down. 

3.2.5. 1  Startup  of  the  Lasers. 

1.  Turn  on  chiller  to  the  Verdi  Nd:YV04.  Allow  it  to  reach  18°C. 

2.  Once  the  chiller  is  at  18°C,  turn  on  the  Verdi  and  open  the  aperture. 

3.  Allow  a  2  minute  warmup  time. 

4.  The  Ti:Saph  should  be  adjusted  for  maximum  power  and  mode-locked.  If  mode¬ 
locking  is  not  possible,  reduce  the  humidity  to  below  25%. 

5.  Turn  on  the  OPO  controller. 

6.  Adjust  the  piezo  voltage  to  1700  and  look  for  a  signal. 

7.  Adjust  the  piezo  voltage  to  get  the  signal  to  1.300  /mi. 

8.  Adjust  the  Ti:Saph  to  maximize  the  power  in  the  signal. 

9.  If  there  is  no  signal,  go  to  the  pump  display  screen  and  make  sure  the  Ti:Saph 
is  set  to  832  nm.  Repeat  6-8. 

3. 2. 5. 2  Startup  of  the  PL  Equipment. 

1.  Get  liquid  nitrogen  (LN2)  from  the  cryo-room. 

2.  Cool  the  monochromator  and  the  cold  finger  with  LN2. 

3.  Turn  on  SRS540  chopper  wheel,  SR850  lock-in  amplifier,  Jobin  Yvon  DataS- 
can2,  IG4400  Ionization  Gauge  Controller,  Alcatel  CFV  vacuum  system  and 
Lakeshore  330  Autotuning  Temperature  Controller. 

4.  Start  the  SynerJY  software. 

5.  Initalize  the  software  to  collect  PL. 

6.  Start  collecting  PL. 
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3. 2. 5. 3  Alignment.  The  first  time  the  system  is  aligned. 


1.  Turn  on  the  Inframetrics  IR  viewer  and  allow  to  initialize.  Ensure  the  cap  is  on 
the  lens. 

2.  Starting  with  either  the  signal  or  idler  path,  use  the  night  vision  goggles  and  IR 
card  to  see  where  the  beam  is. 

3.  Adjust  each  mirror’s  tip  and  tilt  to  steer  the  beam  to  the  next  mirror. 

4.  When  the  beam  is  too  faint  to  see  with  the  night  vision  goggles  due  to  divergence, 
place  a  beam  block  in  the  path. 

5.  Allow  the  beam  block  to  heat  up  and  observe  with  the  IR  viewer. 

6.  Make  small  changes  to  the  previous  mirror’s  tip  and  tilt  and  wait  to  see  the 
thermal  changes. 

7.  Repeat  6  until  the  beam  is  centered  on  the  next  mirror. 

8.  Continue  using  the  IR  viewer  and  beam  block  to  align  the  remaining  mirrors. 

9.  Repeat  steps  3-8  with  the  other  path. 

3. 2. 5. 4  Daily  alignment. 

1.  Place  the  pinhole  in  the  crystal  mount. 

2.  Use  the  photodiode  to  measure  the  irradiance  of  the  idler  and  PL. 

3.  Adjust  the  focusing  mirror  (for  the  PL)  or  lens  (for  the  idler)  to  maximize  signal. 

3. 2. 5. 5  Shutdown. 

1.  Shutdown  is  the  reverse  of  startup. 

3.3  Description  of  the  sample 

This  experiment  will  use  a  p-doped  InAsSb  epitaxially  grown  bulk  semiconduc¬ 
tor  as  the  sample.  It  was  chosen  for  its  large  surface  area,  as  this  allows  the  most 
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leniency  in  alignment.  The  structure  of  the  sample  is  shown  in  Figure  3.2.  The  41- 
period  superlattice  is  used  for  lattice  matching  from  the  p-55%AlGaAsSb  and  the 
p-InAsSb.  The  first  period  has  a  layer  of  123  A  for  the  AlGaAsSb  and  3  A  for  the 
InAsSb.  For  each  period,  the  AlGaAsSb  thickness  decreases  by  3  A  and  the  InAsSb 
thickness  increases  by  3  A.  These  periods  result  in  a  final  layer  of  3  A  of  AlGaAsSb 
and  123  A  of  InAsSb. 


0.5/jmp- InAsSb  (l . 5E 1 6  @77K ) 

=  41  Period  SL 


1( )( )  A  p  -  55%  AlGaAsSb  Field  Stop 
0.5  /an  u  -  55%AlGaAsS b 
0.5  //mn  -GaSb 
n  -  GaSb 


u  -  55%AlGaAsSb 

1 

1  Period 

p- InAsSb 

Figure  3.2.  An  illustration  of  the  structure  of  p-InAsSb  sample 


3.4  Summary 

This  section  contains  a  detailed  look  at  the  equipment  used  in  TRPL.  It  also 
includes  an  enumerated  list  of  procedures  for  starting  up  the  lasers,  starting  the 
equipment,  aligning  the  optics  and  shutting  down  the  system. 
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IV.  Results 


There  were  two  additional  goals  beyond  supplying  documentation  for  the  setup  of 
TRPL  and  operating  procedures.  These  goals  where  to  find  ZPL  and  observe  PL 
from  a  sample  semiconductor.  Since  the  signal  and  idler  beams  were  low  power 
and  difficult  to  align,  ZPL  was  challenging  to  find.  In  addition,  collecting  PL  and 
directing  it  to  the  monochromator  proved  difficult  with  the  low  intensity  PL.  ZPL 
was  attempted  with  a  LiIC>3  crystal  cut  at  20.3°,  and  the  idler  and  signal  beams  at 
ordinary  polarizations  to  the  crystal. 

4.1  ZPL 

One  goal  of  this  experiment  was  to  obtain  ZPL.  This  was  not  obtainable  because 
the  upconversion  efficiency  kept  the  power  of  the  upconverted  beam  below  detection 
with  the  photodiode.  To  obtain  ZPL,  a  PMT  would  need  to  be  used.  Trying  to  obtain 
ZPL  with  the  signal  and  idler  from  the  OPO  consumed  vast  quantities  of  time.  One 
idea  to  improve  the  efficiency  is  to  try  to  find  a  nonlinear  crystal  with  a  higher  de//  as 
the  SFG  efficiency  is  proportional  to  d 2e^.  [23]  This  would  require  close  examination 
of  each  nonlinear  crystal  possible  and  would  be  time  intensive. 

4.2  PL 

A  second  goal  of  this  thesis  was  to  obtain  PL  spectrums  from  a  sample  semi¬ 
conductor.  InAsSb  was  selected  and  its  PL  profile  is  shown  in  Figure  4.1.  The  peak 
intensity  was  shown  to  be  at  3.8  fin 1  (at  ~  150K)  for  this  sample  of  InAsSb.  Roberts 
found  the  peak  to  be  at  3.28  fin 1  (at  80K)  as  shown  in  Figure  4.2.  For  TRPL,  mixing 
for  either  of  these  wavelengths  and  the  1.300  /mi  signal  beam  can  still  be  accomplished 
with  a  LiIC>3  crystal. 

4-3  Varshni  Equation 

To  compare  the  PL  values  of  3.28  and  3.8  /mi  at  two  different  temperatures 
needs  to  take  the  Varshni  equation  into  consideration.  The  Varshni  equation  relates 
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Figure  4.1.  PL  for  InAsSb  from  2000  to  4500  nm.  T  =  150K 
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Sample  106-166,  80K 
peak  at  3280  nm 
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Figure  4.2.  PL  for  InAsSb  from  2400  to  5500  nm  from  Roberts.  T  =  80K 


the  band  gap  energy  to  temperature,  and  is  defined  as 

E„(T)  =  Es(0I<)  - 


(4.1) 


where  Eg  is  the  energy  of  the  band  gap,  T  is  the  temperature  and  a  and  /3  are 
constants  determined  by  the  material.  For  InAs,  Eg(0K )  is  0.415  eV,  a  is  2.76  *  10-4 
^  and  f3  is  83K.  This  gives  units  of  eV  to  Eg{T).  The  corresponding  wavelength  is 
then 


A  = 


h  ■  c 

AT) 


(4.2) 


For  T  =  80K,  A  =  3.07  /irn.  For  T  =  150K,  A  =  3.19  /irn.  However,  the  addition 
of  Sb  would  shift  the  energies.  Taking  3.28  /irn  as  the  correct  value  for  80K,  this 
gives  a  Eg(80K)  of  0.404  eV.  This  is  a  shift  of  -0.026  eV  from  the  value  for  InAs. 
Applying  the  same  shift  to  Eg(150)  yields  0.324  eV.  This  energy  corresponds  to  a  PL 
wavelength  of  3.83  /irn. 


4  •  4  Summary 

PL  from  an  InAsSb  MWIR  sample  at  150K  was  found  to  be  located  at  3.8  /mi. 
This  is  in  agreement  with  Roberts  measurement  at  80K.  ZPL  was  not  found  due  to 
the  time  required  to  align  the  optics. 
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V.  Conclusions 


The  objective  of  this  research  was  to  improve  on  the  technique  for  TRPL.  This  chapter 
summarizes  the  efforts  and  recommends  future  work. 

5. 1  Summary 

TRPL  with  the  signal  and  idler  outputs  from  an  OPO  is  a  difficult  and  time 
consuming  technique  to  measure  carrier  dynamics  of  semiconductors.  Without  a  well 
documented  methodology,  this  work  becomes  almost  impossible  to  accomplish  in  a 
timely  fashion.  Documentation  and  techniques  used  during  previous  efforts  have  been 
changed  by  the  introduction  of  an  OPO.  In  addition,  upconversion  was  attempted  with 
the  1.300-/mi  signal  and  2.311-/un  idler  beams  of  the  OPO  with  a  LTO3  crystal  cut 
at  20.45°  oriented  at  +0.33°. 

PL  data  was  collected  for  a  sample  InAsSb  semiconductor  showing  the  peak 
luminescence  was  at  3.8  //m.  This  wavelength  can  be  mixed  with  the  sampling  beam 
within  the  LiI03- 

5.2  General  Recommendations 

Future  work  would  be  to  continue  the  efforts  at  improving  the  methodology  of 
TRPL  or  to  use  the  available  documentation  to  investigate  the  carrier  dynamics  of 
various  semiconductors.  This  effort  focused  on  the  procedures  for  the  setup  of  TRPL 
whereas  the  continued  improvements  would  be  centered  on  the  collection  and  analysis 
of  data  from  the  luminescence. 

As  stated  in  Chapter  2,  the  use  of  SRH,  radiative  and  Auger  coefficients  can  be 
improved  on  by  expanding  the  rate  equation  and  keeping  the  higher  order  terms  for 
MQW  semiconductors.  For  bulk  semiconductors,  the  use  of  the  SRH,  radiative  and 
Auger  coefficients  is  sufficient  for  realistic  values. 
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Appendix  A.  Software 

This  appendix  describes  the  selection  of  crystals  for  use  in  upconversion.  There  are 
two  steps  involved  in  selection:  use  of  SNLO  and  Mathematica®. 


A.l  SNLO 

SNLO  nonlinear  optics  code  is  a  Windows®-based  program  to  assist  in  selecting 
a  nonlinear  crystal  among  other  features.  One  feature  of  SNLO  is  the  ability  to 
choose  incident  wavelengths  and  available  non  linear  crystals.  It  also  computes  phase 
matching  angles  based  on  polarization  and  the  angle  of  the  crystal  face  to  the  incident 
beams. 

SNLO  is  made  available  by  A.V  Smith  at  Sandia  National  Laboratories  (SNL). 
It  is  no  longer  maintained  by  SNL  but  can  still  be  downloaded  from  the  internet.  [2] 

Table  A.l  shows  a  sample  of  the  SNLO  output  for  the  “Qmix”  module  using 
input  parameters  of  a  LTO3  crystal  and  1.300-/xm  and  2.311-/rm  wavelengths. 


Table  A.l.  Sample  Output  Table  from  SNLO 


2311.0  (0)  +  1300.0  (0)  =  832.0  (e) 

Walkoff  [mrad] 

=  0.00 

0.00 

54.98 

Phase  Velocities  c/ 

=  1.838 

1.852 

1.847 

Group  Velocities  c/ 

=  1.868 

1.875 

1.889 

GrpDelDisp  (fs2/mm) 

=  -46.5 

92.0 

179.0 

At  Theta 

=  20.3 

deg 

deff 

=  1.61E0 

prn/V 

Field  gain  coeff. 

=  6.39E-5 

/sqrt.Watt 

Crystal  ang.  tol. 

=  0.82 

mrad-cm 

Temperature  range 

=  264.25 

K-cm 

Mix  accpt  ang 

=  2.29 

1.28 

mrad-cm 

Mix  accpt  bw 

=  48.15 

74.16 

cm-l-cm 

A. 2  Mathematica 

Once  the  crystal  had  been  selected  via  SNLO,  it  has  to  be  placed  at  an  ori¬ 
entation  that  could  yield  upconversion.  The  use  of  a  numerical  software  package  is 
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necessary  to  determine  the  orientation  of  the  crystal.  A  Mathematica®  notebook  was 
prepared  in  a  previous  dissertation  prospectus  by  LtCol  Roberts  [20]. 

The  calculations  below  describe  the  use  of  this  LiI03  for  an  upconversion  ex¬ 
periment,  where  a  pulsed  beam  at  1300  nm  wavelength  (the  "NL  pump"  beam)  is 
used  to  time-resolve  micl-IR  (2-5  /irn)  photoluminescence  from  a  semiconductor.  The 
semiconductor  is  pumped  by  a  pulsed  1300  nm  beam  (the  "sample  pump"  beam)  that 
has  the  same  temporal  characteristics  (PRF  and  PW)  as  the  NL  pump  beam. 

Initialization 

First,  define  alternate  trig  functions  that  use  degrees  instead  of  radians. 

dSin[x_]  :=Sin[Pix/180] 
dCos  [x_]  :=Cos  [Pix  / 180] 
dTan[x_]  :=Tan[Pix /ISO] 
dArcSin[x_]  :=180ArcSin[x]/Pi 
dArcCos  [x_] : = 180  ArcCos  [x]  /Pi 
dArcTan  [x_] : = 1 80  ArcTan  [x]  /Pi 
e  =  0.00000000001; 

The  index  ellipsoid  for  LiI03  (from  Kato,  same  as  used  by  SNLO)  is: 

ny[A_]:=Sqrt [3.415716  +  (0.047031/(AA2  -  0.035306))  -  0.008801AA2]; 
nz[A_]:=Sqrt[2.918692  +  (0.035145/(AA2  -  0.028224))  -  0.003641AA2]; 

The  transmission  window  for  LiI03  is: 

Almin  =  0.300; 

Almax  =  6.000; 

Sum  Frequency  Generation 

The  wavelengths  we’ll  look  at  are  (in  /irn) : 
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( *  A 1  isthewavelengthtobeupcon  verted  ( PL ,  idler,  orsignal)*) 

A2  =  1.300;  (*  the  NL  pump  wavelength  *) 

A3[A1_]:=1/((1/A1)  4-  (1/A2))(*  the  upconverted  wavelength  *) 

Each  of  the  three  beams  sees  a  different  optical  index: 

nl[Al_]:=ny[Al];  (*or  dinary  beam,  y  —  polarized*) 
n2  =  ny[A2];  (*ordinarybeam.  y  —  polarized*) 

n3[Al_,  0_]:=ny[A3[Al]]nz[A3[Al]]/Sqrt[(nz[A3[Al]]dCos[0])A2  +  (ny  [A3  [  A 1]  ]  dSin  [<9] ) A  2] ; 
(*extraordinarybeam,  xz  —  polarized*) 

where  6  is  the  angle  from  the  z  axis  to  beam  3. 

Next,  set  up  the  phase  matching  requirement: 

kl[Al_]:=2Pinl[Al]/Al; 
k2  =  2Pin2/A2; 
k3min[Al_]:=Abs[kl[Al]  —  k2] 
k3max[Al_]:=kl[Al]  +  k2 
k3[Al_,  0_]:=2Pin3[Al,  0]/A3[Al] 

For  example,  consider  upconversion  of  the  2311  nm  excitation  beam,  scattered 
off  the  sample  (for  ZPL): 

kl[2.311] 

k2 

k3min[2.311] 

k3max[2.311] 

Plot[{k3min[2.311],  k3max[2.311], k3[2.311, 0]},  {9, 0, 180},  AxesLabel  {"0",  "k"}]; 

4.99669 

8.95029 

3.9536 

13.947 
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k 


12; 

10 
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6 

'  25  50  75  100  125  150  175  6 

k3  needs  to  be  between  k3min  =  |kl-k2|  and  k3max  =  kl+k2  for  phase  matching 
to  be  possible.  From  the  plot  above,  we  see  that  this  is  only  true  over  a  limited  range 
of  9: 

0min [A  1  _] : =Replace [x ,  FindRoot[k3[Al,x]==k3max[Al],  {x,  45}]] 
0max[Al_]:=Replace[x,  FindRoot[k3[Al,  x]  ==  k3max[Al],  {x,  135}]] 

Plot[k3max[x]  —  k3[x,  90],  {x, 0.2, 10},  AxesLabel  —*  {"0",  "k3max-k3"}]; 

k3max-k3 
1.5 

1 

0.5 

-0.5 

The  above  plot  shows  the  range  of  PL  wavelengths  for  which  phase-matching 
can  be  achieved,  assuming  all  input  beams  can  be  coupled  into  the  crystal  and  that 
the  crystal  will  transmit  at  the  wavelengths  involved.  Where  the  function  plotted  is 
positive,  phase-matching  is  possible.  The  lower  limit  is: 

Replace[x,  FindRoot[k3max[x]  —  k3[x,  90]  ==  0,  {x,  0.4}]] 

0.285172 
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which  is  outside  the  transmission  window  defined  above,  so  the  lower  Amin  needs 
to  be  reset  to  the  initial  valueis  new  value.  The  upper  limit  is  obviously  outside  the 
transmission  window,  and  so  is  unchanged. 


Continuing  the  ZPL  example: 

0min[2.311] 

0max[2.311] 

20.2703 

159.73 

Note  that  dmin  and  dmax  correspond  to  the  angles  where  k3  is  exactly  equal 
to  kl+k2.  This  is  the  collinear  case,  and  onr  9mm  is  within  rounding  of  SNLO’s 
prediction  of  a  20.3  degree  phase-matching  angle  for  photoluminescence  at  2311  nm. 
As  6  increases  from  dmin,  k3  becomes  less  than  kl+k2,  so  the  angle  between  the  kl 
and  k2  beams  has  to  open  up. 

Now  to  calculate  the  angles  a  and  f3  from  the  first  beam  to  the  second  and  third 
beams,  respectively,  using  the  law  of  cosines. 

a [  A 1  0_] : =d ArcCos [(— kl[Al]A2  -  k2A2  +  k3[Al,  0]A2)/(2kl[Al]k2)] 

/?[A1_,  <9_] : =d ArcCos [(kl[ A  1]A 2  -  k2A2  +  k3[Al,0]A2)/(2kl[Al]k3[Al,0])] 

Plot[o:[2.311,  #],{$,  0min[2.311],0max[2.311]},  AxesLabel  — ►  {"0",  "a"}]; 
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Plot[/?[2. 311,0],  {0, 0min[2.311],0max[2.311]},  AxesLabel  — ►  {"0",  "/?"}]; 
/3[2.311, 90]; 


0 


0 


As  expected,  having  0  =  0min  corresponds  to  the  collinear  scenario,  where 
a— (3— 0.  It  is  also  interesting  to  note  that  the  maximum  spread  between  the  input 
beams  (inside  the  crystal)  that  will  still  allow  phase-matching  is  about  28.5  degrees. 

All  of  the  above  angles  are  internal  to  the  crystal.  The  angles  outside  the  crystal 
need  to  be  investigated,  since  that’s  what  can  be  directly  measured.  For  each  value 
of  0,  there  is  one  and  only  one  pair  ( a}/3 ).  Using  Snell’s  law: 

7l[Al_,  0_,  ^>_]:=dArcSin[nl[Al]dSin[0  —  ip  +  /3[A1, 0]]] 

72[A1_,  0_,  ip-}  :=dArcSin[n2dSin[0  —  ip  +  /5[A1, 0]  —  o:[Al,  0]]] 

73[A1_,  0_,  -0_]:=dArcSin[n3[AU  0]dSin[0  —  U]] 

The  above  equations  have  a  new  variable  ip  to  represent  the  cut  angle  of  the 
crystal.  It  is  the  angle  from  the  z  axis  of  the  crystal  to  the  normal  to  its  surface.  Since 


36 


there  are  limits  on  9,  there  may  be  limits  on  ifj  also.  In  particular,  all  the  external 
angles  need  to  be  in  [-90,90]  so  that  they  will  couple  into/out  of  the  crystal. 

The  upper  limit  on  ifj  is  the  lowest  value  that  makes  either  7I  or  y2  go  to  -90 
degrees;  the  lower  limit  makes  y2  or  7I  go  to  90  degrees.  The  limits  on  ip  are: 

^»min[Al_,  #_]:=#  +  (3[\1,9]  —  dArcSin[l/nl[Al]] 

^»max[Al_,  0_]:=0  —  o:[Al,  6]  +  6]  +  dArcSin[l/n2] 

Plot[{b’rnin[2.311, 9\.  b’rnax[2.311. 0]},  {9,  0min[2.311],  0max[2.311]},  AxesLabel  — >  {"0", 

*l> 

150 

100 

50 

The  above  graph  can  be  interpreted  as  follows.  A  horizontal  line  from  the 
crystal  cut  angle  on  the  vertical  axis  will  intersect  the  two  curves  at  the  minimum  and 
maximum  values  of  9  for  which  all  beams  can  be  coupled  into  and  out  of  the  crystal, 
and  that  allow  phase-matching  for  some  (still  undetermined)  angle  between  the  input 
beams.  Alternately,  for  a  given  desired  phase-matching  angle  on  the  horizontal  axis, 
a  vertical  line  can  be  drawn  that  will  intersect  the  two  curves  at  the  minimum  and 
maximum  values  of  the  crystal  cut  angle  that  will  allow  that  particular  configuration 
to  be  achieved.  Yet  another  interpretation:  All  beams  can  be  coupled  into  and  out  of 
the  crystal,  and  phase-matching  achieved,  for  any  point  (9,i)j)  inside  the  area  bounded 
by  the  two  curves  and  two  vertical  lines  at  dmin  and  dmax. 

Now  to  calculate  the  orientations  of  the  beams  outside  the  crystal.  In  addition 
to  the  7’s  defined  above,  5  is  defined  as  the  spread  between  the  input  beams  for  the 
positive  and  negative  solutions. 

8 [A1-,  0-,  :=7l [A1 72[A1, 9 ,  V>] 
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Using  the  above  equation,  it  is  now  possible  to  find  the  crystal  orientation 
required  to  achieve  phase-matching  for  a  given  wavelength  input  beam  spread,  <5,  and 
crystal  cut  angle,  ip,  by  forcing  the  <5  to  match  the  experimental  setup  for  all  input 
parameters  except  9,  and  then  solving  for  9.  The  solution  can  then  be  used  to  find  y2, 
which  is  the  angle  between  the  normal  to  the  crystal  surface  and  the  direction  of  the 
incoming  NL  pump  beam.  Since  the  direction  and  wavelength  of  the  NL  pump  beam 
are  fixed,  it  makes  a  good  reference  point  for  the  orientation  of  the  crystal.  First  it  is 
necessary  to  "extend"  the  above  definition  of  5  so  that  the  function  will  produce  an 
output  (and  not  an  error)  for  input  values  where  5  is  not  defined  (i.e.  where  phase¬ 
matching  is  not  possible).  Since  6  will  always  be  positive,  extend  the  above  definitions 
by  making  <5=-l  when  it  would  otherwise  not  be  defined,  to  eliminate  ambiguity: 

<fe[Al_,0_,  U-]:=If[(U  >  'i/unin [ A  <  Umax[Al,  0]),5[A1,0,U]>  —1] 

For  SHG  of  the  NL  pump  beam,  the  interaction  is  collinear,  so  9  has  already 
been  calculated:  it’s  dmin: 

0min[2.311] 

Re[t/>min[2.311, 0min[2.311]]] 

Re  [Umax  [2.311. 0min[2.311]]] 

20.2703 

-12.6945 

52.9544 

Again,  the  phase-matching  angle  agrees  with  SNLO’s  prediction  of  20.3  degrees. 
Since  our  crystal  cut  angle  of  20.45  degrees  is  in  the  range  (-12.7,  52.95),  phase¬ 
matching  will  be  possible  and  will  happen  at  the  crystal  orientation: 

Re[72[0.775, 0min[2.311],  20.45]] 

-0.315256 

The  crystal  mount  will  not  block  the  beam  at  this  angle. 
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For  ZPL,  the  spread  angle  becomes  22.4  degrees  and  the  input  wavelength  is 
that  of  the  excitation  beam: 

Plot[<5e[2.311, 0, 39.0]  -  22.4,  {Mmin[2.311],0max[2.311]},  AxesLabel  -►  {"0",  "<5e"}]; 

5e 


It  looks  like  the  solution  is  around  9= 25  degrees: 

x  =  Replace^,  FindRoot[£e[2.311, 9 , 20.45]  ==  22.4,  {(9, 25}]] 

^»min[2.311,  x\ 

^»max[2.311,  x ] 

25.0735 

-0.216952 

53.4821 

Again,  the  crystal  cut  angle  (20.45)  is  within  the  allowed  range  (-0.2,  53.5).  The 
orientation  for  phase-matching  will  be: 

72[2.311,®,  20.45] 

0.644319 

The  direction  of  the  output  (upconverted)  beam  will  be: 

73[2.311f®,  20.45] 

8.51803 

The  direction  of  the  other  input  beam  will  be: 

7l[2.311,3, 20.45] 
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23.0443 


Here’s  a  plot  of  72  vs  a  range  of  values  of  S,  for  2311  nm  wavelength  and  20.45 
degree  crystal  cut  angle: 

Plot[72[2.311,  Replace^,  FindRoot[£e[2.311, 9 , 2.311]  ==  x,  { 9 , 25}]],  20.45],  {x,  19, 25}, 
AxesLabel  —*■  {"0",  "72"}]; 

y2 

-0.6 

-0.8 

-1.2 

-1.4 
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